Effective-one-body waveforms calibrated to numerical relativity simulations: 
coalescence of non-precessing, spinning, equal-mass black holes 
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We present the first attempt at calibrating the effective-one-body (EOB) model to accurate 
numerical-relativity simulations of spinning, non-precessing black-hole binaries. Aligning the EOB 
and numerical waveforms at low frequency over a time interval of lOOOAf, we first estimate the phase 
and amplitude errors in the numerical waveforms and then minimize the difference between numeri- 
cal and EOB waveforms by calibrating a handful of EOB-adjustable parameters. In the equal-mass, 
spin aligned case, we find that phase and fractional amplitude differences between the numerical and 
EOB (2,2) mode can be reduced to 0.01 radians and 1%, respectively, over the entire inspiral wave- 
forms. In the equal-mass, spin anti-aligned case, these differences can be reduced to 0.13 radians 
and 1% during inspiral and plunge, and to 0.4 radians and 10% during merger and ringdown. The 
waveform agreement is within numerical errors in the spin aligned case while slightly over numerical 
errors in the spin anti-aligned case. Using Enhanced LIGO and Advanced LIGO noise curves, we 
find that the overlap between the EOB and the numerical (2,2) mode, maximized over the initial 
phase and time of arrival, is larger than 0.999 for binaries with total mass 30-200Mq. In addition to 
the leading (2,2) mode, we compare four subleading modes. We find good amplitude and frequency 
agreements between the EOB and numerical modes for both spin configurations considered, except 
for the (3,2) mode in the spin anti-aligned case. We believe that the larger difference in the (3,2) 
mode is due to the lack of knowledge of post-Newtonian spin effects in the higher modes. 

PACS numbers: 04.25.D-, 04.25.dg, 04.25.Nx, 04.30.-w 



I. INTRODUCTION 

Coalescing black-hole binaries are among the most 
promising sources for the current and future laser- 
interferometer gravitational-wave detectors such as 
LIGO/ Virgo [IHf and LISA g. 

In general relativity, black holes are defined only by 
their masses and spins; thus generically, a black-hole 
binary depends on eight parameters (mi, Si, m?, S<z)- 
Hence, when black holes carry spins, it is expected that 
tens of thousands of waveform templates may be needed 
in order to extract the gravitational- wave signal from the 
noise using matched-filtering techniques. Considering the 
high computational cost of running numerical-relativity 
simulations of spinning binary black holes (on the order of 
75, 000 CPU hours for 15 orbits, moderate spins and mild 
mass ratios) and the large binary parameter space, it will 
be impractical for numerical relativity alone to provide 
data analysts with a template bank. The work at the 
interface between analytical and numerical relativity 0- 
[IH has demonstrated the possibility of modeling analyti- 
cally the dynamics and the gravitational- wave emission of 
coalescing nonspinning black holes, thus providing data 
analysts with preliminary analytical template families to 
be used for the searches. The next important step is to 
extend those studies to spinning, precessing black holes. 
The present paper represents the first attempt in this 
direction, although limited to non-precessing waveforms, 



within the effective-one- body (EOB) formalism fTM^ 
of spinning black holes [HHj. Recently, Ref. Hf con- 



structed a template family of spinning, non-precessing 
black-hole binaries using a phenomenological approach, 
where the numerical-relativity waveforms are fitted to 
templates which resemble the post-Newtonian (PN) ex- 
pansion, but in which the coefficients predicted by PN 
theory are replaced by many arbitrary coefficients cali- 
brated to numerical simulations. 

The first EOB Hamiltonian that included spin effects 
was worked out in Ref. . In Ref. 22 1 , the authors used 
the non-spinning EOB Hamiltonian augmented with PN 
spin terms to carry out the first exploratory study of the 
dynamics and gravitational radiation of spinning black- 
hole binaries du ring inspiral, merger and ringdown. Sub- 
sequently, Ref. |2(| extended the EOB Hamiltonian of 
Ref. to include next-to-leading order spin-orbit cou- 
plings. In those descriptions, the effective particle is en- 
dowed not only with a mass fx, but also with a spin er. 
As a consequence, the effective particle interacts with the 
effective Kerr background (having spin SKerr) both via 
a geodesic-type interaction and via an additional spin- 
dependent interaction proportional to its spin cr. The 
EOB Hamiltonian developed in Refs. 0, [2(| (with non- 
spinning PN couplings through 3PN order) clarified sev- 
eral features of spinning two-body dynamics. However, as 
we shall discuss below, it is not straightforward to extend 
this Hamiltonian to include higher-order nonspinning 
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PN couplings, such as the 4PN or 5PN adjustable pa- 
rameters recently calibrated to numerical-relativity sim- 
ulations [HI, EH- Moreover, the EOB Hamiltonian of 
Ref. (20j . based on an ad hoc test-particle limit, does 
not reduce to the Hamiltonian of a spinnin g te st parti- 
cle in Kerr spacetime. More recently, Ref. [23[ derived 
the canonical Hamiltonian of a spinning test particle in 
a generic curved spacetime at linear order in the particle 
spin. The construction of an im pro ved EOB Hamilto- 
nian based on the results of Ref. [23[ is currently under 
development. Despite the limitations mentioned above, 
the EOB Hamiltonian of Ref. [2(| is an excellent starting 
point for exploring the calibration of numerical-relativity 
waveforms of spinning black holes within the EOB for- 
malism. Thus, we have used it in this first exploratory 
study, augmenting it with a few adjustable parameters 
that we shall calibrate to two numerical-relativity simu- 
lations. For the EOB non-conservative dynamics, we use 
the gravitational- wave energy flux which includes spin ef- 
fects and which has been computed using the factorized 
multipolar waveforms of Refs. pil l25j. 

The two numerical-relativity simulations we shall use 
describe the evolution of equal-mass, equal-spin, non- 
precessing black-hole binaries. They are produced by the 
pseudospectral code SpEC of the Caltech-Cornell-CITA 
collaboration. In these two configurations, the spins are 
cither aligned ( "up- up" , or UU) or anti-aligned ( "down- 
down" , or DD) with the orbital angular momentum, and 
have dimensionless magnitude Xi = Xi — 0.43655 for 
the UU configuration, and 0.43757 for the DD configura- 
tion. The UU simulation lasts for about 28 gravitational- 
wave cycles or until t = 3250M, and stops about three 
gravitational-wave cycles before merger. The DD simula- 
tion lasts for about 22 gravitational-wave cycles or until 
t = 2500M, and contains the full inspiral, merger, and 
ringdown waveform. Detailed information on the numer- 
ical simulation of the DD configuration can be found in 
Ref. 

For zero spin, the EOB-model considered here agrees 
with the waveform of the equal-mass non-spinning bi- 
nary black hole [53] to a similar degree as the model 
constructed in our earlier work [l5j |. It differs from the 
model presented in [l5| by its modeling of the energy flux 
of gravitational-wave radiation. 

This paper is organized as follows. In Sec. [H] we de- 
scribe the spin EOB model adopted in this paper, includ- 
ing its adjustable parameters. In Sec. IIII1 we calibrate 
the spin EOB model to the numerical simulations, and 
discuss the impact of our results on data analysis. Fi- 
nally, Sec. II VI summarizes our main conclusions. 



II. EOB MODEL FOR SPINNING BLACK-HOLE 
BINARIES 



Xi, Pi, and Si to denote the mass, the position vector, 
the linear momentum vector, and the spin vector of the 
i-th body. We work in the center-of-mass frame defined 
by Pi + P2 = 0. The two body system is described by 
the relative position R = X = X\ — X2 and the relative 
linear momentum P = P\ — —P^. For convenience, we 
define reduced variables 



R P 

M /1 

where M = m\ + 1112 and \i = m\m2/(jrii + 7712). 



(1) 



A. EOB conservative dynamics 

Following Refs. [lil [2(|, we assume that the effective 
particle in the EOB description is endowed not only with 
a mass \i, but also with a spin er. As a consequence, 
the effective particle interacts with the effective Kerr 
background (having spin S^err and mass M) both via 
a geodesic-type interaction and via an additional spin- 
dependent interaction proportional to its spin er. We 
define the Kerr-like parameter a as a = SKerr/M, where 
Skerr denotes the modulus of the deformed-Kerr spin vec- 
tor Skcit- Following Ref. [IH, we write the effective Kerr 
contravariant metric components in a fixed Cartesian-like 
coordinate system. This is done by introducing 
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and writing the contravariant metric components as 
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In this section, we describe the spin EOB model 
adopted in our study and its adjustable parameters. 
Henceforth we use natural units G = c = 1. We use m^, 



1 We denote with the operation of taking the (n, m)-Pade ap- 
proximant. 
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The Taylor approximants to the coefficients A(R) and 
D(R) can be written as 

i=0 

The functions Ak{r) and -Dfe(r) all depend on the sym- 
metric mass ratio v = \xjM through the ^-dependent 
coefficients a,i(y) and di(v). These coefficients are cur- 
rently known through 3PN order (i.e. up to k = 4) and 
can be read off from Eqs. (47) and (48) in Ref. [13| . It 
is worth noticing that although through 3PN order the 
Pade approximant to the function A t (R) of Eq. (|5a|) does 
not pose any problem [2(| , when including 4PN and 5PN 
order coefficients, the Pade approximant develops poles 
for several spin values a. In particular, poles are present 
at large separation when a > 0.7M and the 4PN and 
5PN order coefficient as and are included. 2 Those 
poles could be regularized by adding in Ak(r) higher- 
order spin terms a 2 a^{v)/r 5 , a 2 04(1/) j V 6 and choosing 
for the coefficients 03(1/4) and 04(1/4) negative large 
values (~ —100). Since in this first exploratory study 
we investigate only numerical simulations of moderate 
spins, we do not include any regularization of the poles, 



and consider only the 4PN order coefficient 05. In the 
nonspinning case [t| [l(| EH, EH , the coefficient 05 plays 
an important role in improving the agreement between 
the EOB and numerical waveforms. Here, we choose for 
05 the value obtained by taking the nonspinning limit 
of the spin EOB model and calibrating it to the equal- 
mass black- hole waveform of [27| . following [Tsj . In this 
way, we obtain as(l/4) = 1.775; thus, 05 is no longer an 
adjustable parameter in the spin EOB model. 

In Eq. ((5a)) . we choose m = 1 and n = 4 so that 
A t (R)/R 2 in the limit of a — > reduces to the nonspin- 
ning A(R) used in Refs. [1, E3, EH EH, and we choose 
the same 3PN D{R) function used in those references. 
Therefore, in the spin EOB model, we have 

MR) Num(A t ) 
R 2 Den(A t ) ' [ ' 
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with 
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where x = a/M and a 4 (z/) = (94/3 - 41/32tt 2 ) 1/. Mak- 
ing use of Eqs. ^ and (@|, we can derive 
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The EOB effective Hamiltonian reads [19l \2t 
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L er 



(13) 



where Qi(Pi) is a quartic- momentum term at 3PN order 
independent of spins [HI and L = R x P is the orbital 
angular momentum. In this paper, as a first attempt, we 
use the same spin coupling for the spin a suggested in 
Ref. [20[ , even though it does not reduce to the spinning 
test-particle limit |23| at PN orders higher than 2.5PN. 

In order for H e g to match the PN-expanded spin-orbit 
Hamiltonian through 2.5PN order, we need to require 
that the sum of the spin-orbit couplings of H e s Kerr and 
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-Heffpart gives 

[He&Keri + #cff part] so — ^ L ' (^S^ + 2^ ^ 
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and where the two effective gyro-gravitomagnetic ratios 
gf and gf, read 
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Here a(^) and b(y) are two gauge parameters related to 
the freedom of applying a canonical transformation in- 
volving spin variables. If we knew the exact Hamilto- 
nian, the choice of these parameters should not affect 
the physics of the EOB model. However, since we start 
with an approximate Hamiltonian that reproduces the 
spin-orbit couplings only through 2.5PN order, we expect 
the EOB model to depend on the choice of a(v) and b(v). 
Considering the structure of the gyro-gravitomagnetic ra- 
tios, such dependence should start at 3.5PN order as a 
spin-orbit coupling term. Because of this dependence, 
a{v) and b{v) can be used as adjustable parameters. 

Moreover, in order for H c s to match the PN-expanded 
spin-spin Hamiltonian through 2PN order, the simplest 
choice is to require that the Kerr spin [20j 



"Kerr — S -\- S* . 

As a consequence, Eq. (fT4|) implies 



a = hgf-2)S + hgf,-2)S* 



(17) 



(18) 



To include higher-order spin-spin contributions in the 
EOB effective Hamiltonian, we introduce a 3PN spin- 
spin term whose coefficient a|g N is currently unknown 
and can be used as an adjustable parameter 

H c s(R, P, Si, S 2 ) = H c g KciI (R,P,S 

Kerr ) 

+ H eSpait (R,P,<T) 

+ a 3PN y *Ke^ (ig) 



Finally, the EOB Hamiltonian is 



H, 



real 



Mc 




(20) 



In summary, in this first exploratory study, we choose 
to employ only two adjustable parameters 3 : b(v) which 
introduces a spin-orbit term at 3.5PN order, and a|g N 
which introduces a 3PN spin-spin term. As we shall see, 
these two adjustable parameters are sufficient to reduce 
the phase and amplitude differences between EOB and 
numerical waveforms of the UU and DD configurations 
to (almost) the numerical error. The remaining flexibil- 
ity of the spin EOB model can be exploited in the future 
when numerical relativity simulations of other spin con- 
figurations will become available. Thus, for the rest of 
the paper, we set a(v) in Eq. (|16a[) to zero. 

Within the Hamiltonian approach, radiation-reaction 
effects can be incorporated into the dynamics in the fol- 
lowing way [13, HI: 



dX l . 9-Hreal 

{A , -HrealJ- 



dt 

dP, 



(21) 



= {Pi, H ieal } + f 1 = + ^ ■ ( 22 ) 

Here, Ft denotes the non-conservative force, which is 
added to the evolution equation of the (relative) momen- 
tum to take into account radiation-reaction effects. Fol- 
lowing Ref. [iH, we use 4 



Fi 



1 dE 



Q\L\ dt 



P. 



(23) 



where fl is the orbital frequency and L is the orbital 
angular momentum. The gravitational-wave energy flux 
dE I dt is obtained by summing over the gravitational- 
wave modes (l,m) as 



dE 1 V- k 

dF" 16^ 1^ \ h t™ 



1=2 i 



(24) 



which reduces to the following expression for circular 
equatorial orbits in the adiabatic approximation: 



dE 

~dt 16tt 



8 I 



(25) 



1=2 i 



where Cl is the reduced orbital frequency Ct = MCI. We 
shall define the EOB waveforms htm in Sec. Ill CI The 



3 We find that b(u) is strongly degenerate with a(u). 

4 We notice that this choice of the radiation-reaction force intro- 
duces a radial component of the force R ■ F oc R ■ P = R Pr . 
In the nonspinning EOB models, this component is usually ig- 
nored Q3, 
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equations of motion for the spins are simply obtained 
through the equations 

^i = {Si,M = ^x Sl , (26) 
at obi 

±S 2 = {S 2 , ff roal } = x S 2 . (27) 

at ob2 

In the nonspinning case, it is useful [ljl [H| to re- 
place the radial momentum Pr with Pr. , the conjugate 
momentum of the EOB tortoise radial coordinate R*: 
dR* /dR = \f~D I A. This replacement improves the nu- 
merical stability of the EOB equations of motion because 
Pr diverges when approaching the zero of A(r) (the EOB 
event horizon) but Pr* does not. Therefore, in the spin- 
ning EOB Hamiltonian, we similarly choose to use the 
conjugate momentum to the tortoise radial coordinate of 
the i/-deformed Kerr geometry: 

dR*_ = R 2 + a 2 _J_ 

dR VATA^ UR)' 

In the limit a — > 0, Eq. (|28[1 reduces to the nonspin- 
ning EOB tortoise coordinate defined above. In the limit 
v — > 0, Eq. ((28)) reduces to the tortoise coordinate of 
the Kerr geometry: dR* /dR = (R 2 + a 2 )/ A. Since the 
EOB Hamiltonian and Hamilton equations are written in 
Cartesian coordinates, some algebra is needed to rewrite 
them to include this transform of the radial coordinate. 
In Appendix El we write down explicitly the transformed 
EOB Hamiltonian and Hamilton equations in Cartesian 
coordinates. In particular, Eqs. (f2Tj) and ([22| should be 
replaced by Eqs. (|A"4|) and (|A"5j) . 

Initial conditions for the Hamilton equations are con- 
structed following the prescription of Ref. (22j, which 
provided post-circular initial data for quasi-spherical or- 
bits when neglecting spin-spin and next-to-leading order 
spin-orbit effects. Note that exact circular orbits cease 
to exist in the conservative dynamics when spin-spin and 
next-to-leading order spin-orbit effects are present, ex- 
cept for special configurations in which the spins are 
aligned or antialigned with the orbital angular momen- 
tum. We start each evolution at a large initial separation 
of 50M. The EOB trajectory is sufficiently circularized 
when reaching a separation of ~ 16M, where numerical 
waveforms start. In this way, we remove the residual ec- 
centricity in the EOB trajectory due to imperfect initial 
conditions, while physical eccentricity due to spin effects 
is preserved. 

As a final remark, the spin variable in the EOB model 
is the constant spin variable, i.e., its magnitude does not 
change during precession [29|. We identify it with the 
spin variable in the numerical simulation, which also re- 
mains constant during the evolution (26j . 



B. Characteristics of EOB orbits for spinning, 
non-precessing black holes 

Here we investigate certain properties of the spin EOB 
Hamiltonian that are crucial when building the complete 
EOB model. Specifically, we check the existence and be- 
havior of the innermost stable circular orbit (ISCO), the 
light ring (photon orbit) and the maximum of the EOB 
orbital frequency. This study is restricted to circular 
equatorial orbits in the spin aligned or anti-aligned cases. 
For convenience, we consider the EOB Hamiltonian writ- 
ten in spherical coordinates; we fix 9 = n/2 and set the 
conjugate momentum Pg = 0. 

The ISCO position is obtained by solving dH(R, Pr. = 
O,P )/<9P = 0,d 2 H{R 7 P R , = O,P )/<9.R 2 = where 
Pr* and P^ are conjugate momenta of the tortoise radial 
coordinate R* and the orbital phase <j>, respectively. In 
the nonspinning limit, we find the following ^-correction 
of the ISCO frequency relative to the Schwarzschild value 

^isco = 6~ 3/2 x 

[l + 0.9837^ + 1.2543j/ 2 + 5.018v s + 0{is 4 )] , 

(29) 

where 0(v 4 ) terms contribute less than 1% to the correc- 
tion. In the test-particle limit, the coefficient of the linear 
i/-correction term, 0.9837, should be compared to the re- 
cently available self-force result [3(| (transformed to the 
gauge condition and mass convention used in the EOB 
formalism by Ref. [3l|) of 1.2513. The relative differ- 
ence of 21% is due to the fact that our nonspinning EOB 
Hamiltonian, although calibrated to equal-mass numer- 
ical simulations, does not capture all the ^-dependence 
correctly at 4PN order. 5 The improved spin EOB Hamil- 
tonian [32| will incorporate consistently the self-force re- 
sult (e.g., Ref. [3l|) and can be better constrained by new 
numerical simulations. 

In the spin aligned or anti-aligned case, we find that 
the ISCO exists for all spin magnitudes. However, in 
the spin aligned case, when a > Q.8M, the ISCO radius 
(frequency) starts to increase (decrease) with increasing 
a. This is contrary to the monotonic dependence of the 
ISCO radius (frequency) on the spin magnitude in the 
test-particle limit. This unusual behaviour will be over- 
come by the improved spin EOB Hamiltonian of Ref. {32| . 
Nevertheless, since this problem occurs only at extreme 
spin magnitudes and here we have numerical waveforms 
of moderate spins (\a\ < 0.5M), we choose to use this 
spin EOB Hamiltonian in the current calibration. 



5 We notice that if we used the 4PN and 5PN coefficients, as and 
(26 , suggested in Ref. |3ll , we would obtain poles in the function 
At(R) for \a\ > 0.75M. Moreover, if we adopted the values of as 
and aa , suggested in Ref. [3ll | for the spin configurations analyzed 
in this paper, for which there are no poles in At(R), wc would 
obtain phase disagreements on the same order of the ones we 
have found. 
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The light ring is the unstable circular orbit of 
a massless particle (such photon) and can be 

computed from the deformed EOB metric or from 
H e fiKerr(R, P, Skch)- As in the nonspinning case, we 
do find a light ring with our spin EOB Hamiltonian. 
However, in contrast to the nonspinning case, for several 
values of the spin parameters (including the DD con- 
figuration) our spin EOB Hamiltonian does not yield a 
maximum orbital frequency. It is worth mentioning that 
if we were using only the "Kerr" part of the spin EOB 
Hamiltonian, i.e. we ignore H e g part (R,P,cr), then we 
do find a maximum of the orbital frequency and its value 
is quite close to the light ring position. A more detailed 
study has revealed that the absence of the maximum of 
the orbital frequency for the full spin EOB Hamiltonian 
is due to the spin-orbit coupling term H e ff par t(i2, P, er) 
defined in Eq. (fT3|). which as discussed above does not 
reduce to the test-particle limit prediction at PN orders 
higher than 2.5PN. In the improved spin EOB Hamilto- 
nian [32j], H e ff palt (R,P,cr) will be consistent with the 
test-particle limit prediction at all PN orders linear in 
the particle s pin . Analyses using the improved spin EOB 
Hamiltonian [32J have shown a reasonable agreement be- 
tween the position of the EOB light ring and the maxi- 
mum of the EOB orbital frequency. 

Quite interestingly, when the numerical and EOB 
waveforms are aligned at low frequency, as discussed in 
detail in Sec. IfflAl we find that the EOB light ring is 
reached at time 0.3M before the peak of the numerical 
/122 amplitude. Therefore, a nice property of the nonspin- 
ning EOB model [l4|, [HI holds also in the spinning case, 
i.e. the EOB light ring position is a good approxima- 
tion of the peak position of the numerical /122 amplitude. 
The latter property will be a key ingredient in the EOB 
waveform model, as described later in Sec. Ill Dl 



C. EOB waveform: Inspiral &: Plunge 

Having described the inspiral dynamics, we now turn 
to the gravitational waveforms hg m . The latter can also 
be employed to compute consistently the inspiral dy- 
namics through the radiation- reaction force [l4j. In the 
nonspinning case, Refs. [HI, [H|] have shown that the re- 
summed, factorized PN waveforms proposed in Ref. [25j] 
are in excellent agreement with the numerical waveforms. 
In Ref. [24[ we have generalized the resummed factorized 
waveforms to include spin effects. 

The resummed waveforms are written as the product 
of five factors, 



,(N,e) Me) rp 



(30) 



where e denotes the parity of the multipolar waveform. In 
the circular- orbit case, e is the parity of I + m: e = 7r(l + 
m). These factors are discussed extensively in Ref. [25j . 
Here we simply write down the expressions used in our 
spin EOB model, valid for spins aligned or anti-aligncd 
with the orbital angular momentum. Thus, we restrict 



ourselves to the equatorial plane (6 = ir/2 and po = 0). 
The leading term h^J^ is the Newtonian contribution 



AN,, 



Mv 



n%c i+ ^)v^yt-e,- m ^), (31) 



where 1Z is the distance from the source. The n^L and 
ct+eiv) are functions given in Eqs. (5)— (7) of Ref. [25[. 
The y (0,(f>) are the scalar spherical harmonics. The 
tangential velocity v^, is the non-Kcplerian velocity of a 
spherical orbit defined by = r n Ct where 



-2/3 



M 



dH cS (P R = 0) 



dPs 



-2/3 



(32) 

and -P^.cir is the solution of the spherical orbit condition 
dH c ff(R 7 Pr = Q,P$)/dR = 0. As in the nonspinning 

case, the functions S^g, Ti m , e lSlm and fi m appearing 
in the right hand side of Eq. (|30[) are computed using 
the Keplerian velocity v — Ct 1 / 3 . Moreover, S^} is an 
effective source term that in the test-particle, circular- 
motion limit contains a pole at the EOB light ring. It is 
given in terms of the EOB dynamics as 



S, 



(+) 



A-) 



r x p 



(33) 



Setting S^g to \r x p* \ in Eq. (|33]l is not the only possi- 
ble choice; for example, one may instead choose S l 



(-) 

efi 



to 



be either H eS or J cS — \r x p* + S Kc rr/ (M The ef- 
fect of this choice on the spin EOB model investigated in 
this paper is marginal, since in the equal-mass, equal- 
spin, non-prccessing binary configurations, odd parity 
modes contribute only a tiny fraction of the total en- 
ergy flux (see Sec. IIII Dl for details). Although we choose 
to use the source term defined in Eq. (|3"3"j) . there is no 
evidence indicating that this choice is better or worse 
than others for those binary configurations in which odd 
parity modes are more important. 

The function Ti m in the right hand side of Eq. (f3"0)) 
resums leading logarithms of tail effects, and e lStm is a 
phase correction due to subleading logarithms. Through 
2PN order, there are no tail contributions due to spin 
effects and Tg m and e lSlm do not differ from the non- 
spinning case. Their explicit expressions are given in 
Eqs. (19)-(29) of Ref. j25|. Finally, the functions fi m in 
the right hand side of Eq. (|3"D|) collect the remaining PN 
terms. We computed (24| the spin terms in fe m by Taylor 
expanding the he m in Eq. (|30|) and comparing it to the 
Taylor-expanded ht m calculated in PN theory, including 
the test-particle spin effects through 4PN order. In the 
test-particle limit, we choose S'koit as the spin variable 
of th e sp acetime. Expressions of fe m can be read from 
Ref. m 



For odd parity modes, depending on the choice of the source 
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[15l |25| . we resum all the nonspin- 



in the functional form /JJ^ = (pe m Y 



Following Rcfs. 
ning terms in f tm m mc luu^iuutu ^ TO 
that holds at known PN orders, where /J^ collects the 
nonspinning terms mff m , and pi m can be read from 
Appendix C of Ref. |25|]. The motivation for this p- 
resummation is to reduce the magnitude of the 1PN co- 
efficients in fi m that grow linearly with I (see Sec. IID 
of Ref. [25| ) . Since at leading order we did not find such 
dependence on I in the spinning terms J24[, we do not 
apply the p-resummation to them. 

Furthermore, since we are trying to reproduce effects 
in the numerical simulations that go beyond the quasi- 
circular motion assumption, motivated by the PN ex- 
pansion for generic orbits, we include non-quasicircular 
(NQC) effects in hi m in the form 



hZ p - plunsc ^ NQC h em = h im 

P 2 r , 1 



1 I II 

1 + a t 



h lm Pr 



(rfi) 2 r 

Pr* 1 
{rCl) 2 r 2 



(rf}) 2 r 3 / 2 



(34) 



D. EOB waveform: Merger & Ringdown 

The merger-ringdown waveform in the spin EOB model 
is built in the same way as in the nonspinning EOB 
model. Details on building merger-ringdown waveforms 
can be found in Sec. IIC of Ref. [Hj|. Here we briefly 
summarize the key points. 

In the spin EOB model, the ring-down waveform is a 
linear combination of eight quasinormal modes (QNMs) 
of the final Kerr black hole. Mass and spin of the fi- 
nal black hole are computed from numerical data. In 
particular, for the numerical simulation of the DD con- 
figuration, we use M BK /M = 0.961109 ± 0.000003 and 
»bh/A/bh = 0.54781 ± 0.00001 computed in Ref. 
Frequencies of the QNMs are computed by interpolat- 
ing data from Ref. [3j]. The eight complex coefficients 
of the linear combination are fixed by the hybrid comb 
matching described in Sec. IIC of Ref. [l5j . The 
matching time i^a tch is fixed to be the EOB light-ring 
position. The matching interval Ai^™ tch is an adjustable 
parameter that is fixed by reducing the difference against 
numerical merger-ringdown waveforms. 



III. CALIBRATING THE EOB WAVEFORMS 
TO NUMERICAL RELATIVITY SIMULATIONS 



A similar expression was used in Ref. [l5| except that 
there we used r instead of p r * . For a test-particle plung- 
ing in the Kerr geometry, r goes to zero at the horizon. 
We observe a similar behavior in the EOB ^-deformed 
Kerr geometry. Therefore, in contrast to the nonspin- 
ning case, the evolution of r is not monotonic during the 
inspiral-plungc: r increases during the inspiral, reaches a 
peak, and then starts decreasing during the plunge. By 
replacing r with p r * , we keep the NQC correction terms 
in Eq. Q34p monotonic in time; thus, they can success- 
fully model the monotonically increasing amplitude dif- 
ferences between the quasi-circular EOB and numerical 
waveforms. As in Ref. [l5j . we fix two of the four ad- 
justable parameters a i 22 by requiring that the peaks of 
the numerical and EOB /122 waveforms coincide in both 
time and amplitude, where the peak time of the numer- 
ical /122 waveform is accurately predicted by the EOB 
light ring, as discussed above. The other two a^ 22 param- 
eters are determined by minimizing the overall amplitude 
difference with respect to the numerical waveform as ex- 
plained in detail below. The NQC corrections in Eq. (fM)) 
also depend on spins. However, there is not enough nu- 
merical information in this work (we have only the DD 
configuration) to discriminate between the spinning and 
nonspinning contribution. 



We now calibrate the spin EOB model against the nu- 
merical UU and DD spin configurations. We extract 
both the Newman-Penrose (NP) scalars $4 and the 
strain waveforms hg m from the simulations. The strain 
waveforms are extracted with the Reggc- Wheeler- Zerilli 
(RWZ) formalism (see Appendix of Ref. [lj] for 

details of the numerical implementation used to obtain 
him). We use the RWZ hi m to calibrate the EOB model, 
and use the NP ^f" to check the consistency between the 
two wave-extraction schemes and to estimate the numer- 
ical error associated with the wave extraction schemes. 

We will use the i = 2, to = 2 component of the nu- 
merical waveform for tuning the EOB model. Thus, we 
calibrate in total the following six adjustable EOB pa- 
rameters: b{u), af s N 7 a^ 22 , a^ 22 , a'^ and A^ atch . 



A. Uncertainties in numerical waveforms 

In this section, we compare numerical waveforms com- 
puted at different numerical resolutions and/or using dif- 
ferent extrapolation procedures, or with different wave- 
extraction schemes. Estimates of numerical errors in the 
waveforms will set our standards when calibrating the 
EOB model. 

First, we adopt the same waveform- alignment proce- 
dure used in Ref. [l5| , that is we align waveforms at low 
frequency by minimizing the quantity 



term among -ff e ff , ^eff and j e ff , the corresponding choice of 
should be made among the expressions of fijL , /£l and // 



H(At,A0) 



[&(*)- 02 (t- At)- A(j>} 2 dt, (35) 
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FIG. 1: (color online). Numerical error estimates for the UU 
configuration. We show the phase difference between several 
numerical ^P 22 waveforms aligned using the procedure defined 
by Eq. (|35j>. 
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FIG. 2: (color online). Numerical error estimates for the DD 
configuration. We show the phase difference between several 
numerical ^P 22 waveforms aligned using the procedure defined 
by Eq. (|35). 



N5 (as opposed to N6) is from a simulation with a lower 
numerical resolution and is used to quantify the uncer- 
tainty due to numerical truncation errors. The waveform 
labeled by "h(N6, n=3)" is generated by twice differen- 
tiating the RWZ-extracted "/i(N6, n=3)" waveform, and 
is used to quantify the uncertainty due to the systematic 
difference between extracting the NP scalar and extract- 
ing the strain waveform via the RWZ formalism. 

The noise before t = 500M is due to spurious radiation 
from initial conditions. The features around t « 2I00A/ 
in Fig. [2] are due to a change of gauge in the numerical 
simulation. Extrapolation with n = 2 leads to systematic 
errors in the extrapolated waveform which in turn results 
in a systematic error in At. Therefore, the green dashed 
lines in Figs. [T] and [2] represent a possibly overly conserva- 
tive error estimate. There is a tiny frequency difference 
between the NP and RWZ extracted waveforms, which 
is magnified into a substantial time shift when the wave- 
forms are aligned at low frequency. As a consequence, 
the dot-dashed brown line in Fig. [5] shows a larger phase 
difference which builds up during the late inspiral. It 
provides us with the most conservative error estimate 
for the DD configuration. This is better illustrated in 
Figs. [3] and 01 where we compare the RWZ /122 and NP 
, F 22 waveforms without any time or phase shift. In blue 
solid lines, we show the phase and relative amplitude dif- 
ferences A</>np and AA^p/A between the RWZ /122 wave- 
form differentiated twice with respect to time and "J 22 . 
In red dashed lines, we show the phase and relative am- 
plitude differences A<^rwz an d AA^wz/A between \E' 22 
integrated twice in time and the RWZ h 2 2- In Fig. 2J 
A</>np shows a slope between t = 500M and 1500M. 
When we apply the alignment procedure, this slope is 
removed through a time shift, which is transformed into 
a larger phase difference during late inspiral where the 
wave frequency is large. 



B. Calibrating the EOB adjustable parameters 



over a time shift At and a phase shift A<f>, where <f>i(t) 
and 4>2{t) are the phases of the two waveforms. The range 
of integration (£1,^2) is chosen to be as early as possible 
to maximize the length of the waveform but late enough 
to avoid the contamination from junk radiation present 
in the numerical initial data. The range of integration 
should also be large enough to average over numerical 
noise. We fix t x = 500M and t 2 = 1500M in Eq. (1551) . 

Using this alignment procedure, we estimate the errors 
on the numerical , F 22 . Figs. [TJ and [2] summarize the 
phase errors for numerical . The numerical waveform 
labeled "(N6, n=3)" [or "* 4 (N6, n=3)"] is the reference 
numerical waveform used throughout this paper. Each 
waveform is extracted on a set of spheres at fixed dis- 
tances from the source, and then extrapolated to future 
null infinity; the labels n refer to different orders of this 
extrapolation and are used to quantify the uncertainty in 
the phase due to extrapolation. The waveform labeled by 



Here we adopt the procedure suggested in Ref. [l5l ]. 
augmented with the iterative scheme suggested in 
Ref. [3] when calibrating the adjustable parameters. 

We divide the adjustable parameters into three groups 
and calibrate them in two steps. The first group, EOB- 
dynamics parameters, consists of b(v) and a|s N in the 
EOB Hamiltonian (there is no adjustable parameter in 
the model of the EOB energy flux). These parameters 
determine the inspiral and plunge dynamics of the EOB 
model and affect the merger-ringdown waveform only 
indirectly through the waveform's phase and frequency 
around the matching point. The second group, EOB- 
NQC parameters, consists of a^ m , which enter both 
the EOB dynamics (through the energy flux) and the 
EOB waveform (through the NQC correction). The third 
group, EOB-waveform parameters, consists of Ai^ tch , 
which affect the EOB mcrgcr-ringdown waveform but 
not the EOB inspiral-plunge waveform. All the EOB 
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FIG. 3: Phase and relative amplitude difference between the 
(I, m) — (2, 2) modes of the RWZ waveform hnwz and NP 
scalar ^4 for the UU case. 
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FIG. 4: Phase and relative amplitude difference between the 
(I, m) — (2, 2) modes of the RWZ waveform /irwz and NP 
scalar ^4 for the DD case. The right panel shows an enlarge- 
ment of merger and ringdown, with the dotted vertical lines 
indicating time of maximum of |^4|, and where | M/ | has de- 
cayed to 10% and 1% of the maximal value. (The blue lines 
are smoothed; the grey data in the background represents the 
unsmoothed data.) 



adjustable parameters are calibrated to the numerical 
RWZ h-22- In the first step of calibration, we simulta- 
neously reduce the difference in waveforms against the 
numerical UU and DD configurations by calibrating the 
EOB-dynamics and the EOB-NQC adjustable parame- 
ters. In the second step, using the adjustable parameters 
fixed in the first step, we calibrate the EOB-waveform 
adjustable parameters. 

We adopt the iterative scheme suggested in Ref. [l4[ 



to fix the EOB-dynamics and the EOB-NQC param- 
eters in the first step of calibration. In each it- 
eration, we first minimize the phase difference by 
calibrating the EOB-dynamics adjustable parameters. 
To be specific, we choose to minimize the quantity 
max [0eob(*) - 0NR.(i)] ~ min [0eob(*) - 0nr(<)]- By 
comparing the EOB model to both the UU and DD 
configurations it is possible to calibrate the parameters 
b(v) and a|g N separately. This is because b{y) alters 
a 3.5PN spin-orbit coupling term that depends on the 
spin orientation, so the phases of the UU and DD wave- 
forms change in opposite directions when varying b(v), 
but a|g N alters a 3PN spin-spin coupling term, so the 
phases of the UU and DD waveforms change in the same 
direction when varying &gg N - 

The EOB-NQC adjustable parameters are calibrated 
only to the numerical /122 waveform of the DD config- 
uration, because we did not run the UU case through 
merger and ringdown. We first fix a\ 22 and a^ 22 by re- 
quiring that a local extremum of the EOB hii amplitude 
coincides with the peak of the numerical h^2 in time and 
amplitude (the peak time is predicted by the EOB light 
ring). We expect that in the future, the peak amplitude 
of the numerical /122 will be predicted by numerical rela- 
tivity with high accuracy as an interpolation function on 
the physical parameters. Therefore, a\ 22 and 4 122 can be 
determined without a least-squares fit to the NR wave- 
form, reducing by two the number of parameters to be 
determined by a least-squares fit. The other two NQC 
parameters, aj 22 and a^ 22 , are calibrated to the numeri- 
cal waveform to further reduce the disagreement in am- 
plitude. The NQC parameters will enter the flux through 
the NQC waveform N( 3 c /i22 in the next iteration. They 
are set to zeros initially to start the iteration and they 
usually converge within five iterations. 

In the third step, we calibrate the EOB-waveform ad- 
justable parameter A<^ atch by reducing the difference in 
the DD configuration merger-ringdown /122 waveform. 



C. Comparing the gravitational- wave modes ^22 

Before calibrating the EOB adjustable parameters, we 
investigate the phase difference for the EOB uncalibrated 
waveforms. For the uncalibrated model, we set 05 = 
1.775 and all six of our adjustable parameters to zero. 
We find that during the inspiral, the phase agreement 
between the numerical and spin EOB uncalibrated wave- 
forms is already substantially better than the agreement 
between numerical and Taylor-expanded PN waveforms. 
For the latter, we consider the 3.5PN spin Taylor model 
(T4) of Ref. [HI with amplitude corrections through 
the highest PN order currently known [24], In fact, 
using the uncalibrated spin EOB model and aligning the 
waveforms with the procedure defined by Eq. (|35p . we 
find that the phase differences against the numerical UU 
and DD waveforms, at the end of the simulation and at 
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FIG. 5: (color online). Comparison between the numerical 
and EOB waveform for the UU configuration using b(v) = 
— 1.65 and als N = 1-5- The top panels show the real part 
of the numerical and EOB /122, the bottom panels show am- 
plitude and phase differences between them. The left panels 
show times t = to 2950M, whereas the right panels present 
an enlargement of the later portion of the waveform. 
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FIG. 6: (color online). Comparison between the numerical 
and EOB waveform for the DD configuration using b(v) = 
— 1.65 and a|s N ~ 1-5- The top panels show the real part 
of the numerical and EOB /122, the bottom panels show am- 
plitude and phase differences between them. The left panels 
show times t = to 2300M, and the right panels show times 
t = 2300M to t = 2480A/ on a different vertical scale. 



merger, respectively, are —0.2 and 4.3 rads 7 . Using the 



7 If in the uncalibrated EOB model, we chose to include both 
05 and 0,6, as discussed in Sec. Ill Al and adopted the values 
05 = —15.5 and = 223 (calibrated to equal-mass nonspinning 
numerical waveforms and consistent with the constraint derived 
from self- force results in Ref. [3ll] ). we would find for the phase 
differences —0.3 and 3.5 rads. They are comparable with the 
differences found in the spin EOB model with only as . 



spin Taylor T4 model, the corresponding phase differ- 
ences are 2.0 and —10.0 rads. Therefore, the spin EOB 
model, even uncalibrated, improves the phase agreement 
with numerical waveforms of Taylor-expanded PN mod- 
els by resumming the PN dynamics. 

When calibrating the EOB model, we find that the 
difference in phase and amplitude between the numerical 
and EOB waveforms is minimized when fixing the EOB- 
dynamics parameters b{v) = —1.65 and a|g N = 1-5. 

In Fig. [5l we compare numerical and EOB hii wave- 
forms for the UU configuration. The phase difference and 
relative amplitude difference arc strictly within 0.01 rads 
and 1%, respectively. The systematic error in the EOB 
waveform in the UU configuration is therefore smaller 
than the numerical errors. 

In Fig. [SI we compare numerical and EOB /122 wave- 
forms for the DD configuration. Using b(y) = —1.65 
and a|g N = 1.5 again, we find that the best phase 
and amplitude agreement is obtained when the match- 
ing occurs at an interval of Ai^ atch = 3.5M ended at 
*match = 2402. 0M, which is the EOB light-ring position 
and is 0.3M before the merger, i.e. the peak of the nu- 
merical at t = 2402. 3M. The NQC parameters are 
a\ 22 = -16.1052, 4 22 = -1124.43, aj 22 = 4529.21 and 
a^ 22 = —4587.53. The relative amplitude difference is 
strictly within 1% until 2000M. After 2000M, although 
oscillations due to numerical gauge effects in the RWZ 
/122 waveform are at the level of 2% until the merger, 
the average difference is still less than 1%. After the 
merger, the amplitude difference grows to about —5% 
and starts oscillating with increasing magnitude. The 
latter phenomenon is due to gauge effects in the RWZ /122 
waveform as discussed in Sec. IIII Al and the Appendix of 
Ref. [l~5|. The phase difference is within 0.01 rads until 
about 1800M and grows to —0.28 rads until merger and 
settles to about 0.1 rads before the exponentially decay- 
ing amplitude results in increased errors in the extracted 
gravitational- wave phase. 

In the spin DD configuration, the error in the EOB 
waveform (especially in the phase) is not within the nu- 
merical errors. The phase difference of ~ 0.13 rads at late 
inspiral around t = 2300M can be reduced to within the 
numerical errors of ~ 0.01 rads by calibrating the EOB- 
dynamics adjustable parameters, i.e. b{v) and flgg N - 
However, this leads to an increase of the phase differ- 
ence around the merger. Since we choose to minimize 
the span of the phase difference over the range of the full 
inspiral, merger and ringdown waveform, the phase dif- 
ference at late inspiral is larger than what it could have 
been. The largest phase difference around merger can 
not be removed by calibrating the chosen adjustable pa- 
rameters. Nevertheless, we can substantially reduce the 
phase difference if we allow one of the EOB-dynamics pa- 
rameters b{v) and a|g N to be different in the UU and DD 
cases, or if we add one more spin-independent adjustable 
parameter. For instance, there can be a NQC correction 
factor on the right-hand side of Eq. ([3"4"]) that contributes 
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solely to the phase of the waveform in the form of [28 



/ 4m P ~ plUngC = NQ % m e i6 -*" p^/(^1 (36) 



We can reduce the phase difference at merger to < 0.05 
rads by calibrating this extra EOB-NQC adjustable pa- 
rameter b\ lm . However, until we study a larger number 
of waveforms, we will not over-tune parameters, since the 
main purpose of this preliminary and exploratory study 
on the spin EOB model is to show that by making a very 
simple and minimal choice of adjustable parameters, we 
can achieve a quite fair agreement with the numerical 
simulations. 

We shall emphasize that, despite the small phase dif- 
ference that exceeds the numerical errors in the DD con- 
figuration, the faithfulness of the EOB waveforms with 
the numerical waveforms is very good. Using the noise 
curves of Enhanced LIGO and Advanced LIGO 8 , for both 
the UU and DD configurations, we find that the faithful- 
ness is always better than 0.999 for black- hole binaries 
with a total mass of 30-200M Q . Note that the numerical 
waveforms start roughly at 40Hz for binaries with to- 
tal mass 30M Q and at 10Hz for binaries with total mass 
100M Q . Since the Advanced LIGO noise curve has a 
low frequency cutoff at 10Hz, the numerical waveforms 
arc not long enough to cover the Advanced LIGO sen- 
sitivity band for binaries with total mass smaller than 
IOOA/q. When computing overlaps for these lower mass 
binaries using Advanced LIGO noise curve, we start the 
integration at the starting frequency of the numerical 
waveforms instead of 10 Hz. For the available numerical 
waveforms, the overlaps with EOB waveforms are well 
above the requirement on the accuracy of binary black- 
hole waveforms for detection and measurement purposes 
in gravitational- wave observations [40l ] . 

In Fig. we compare the amplitude and frequency 
of numerical and EOB h^i waveforms together with the 
orbital frequency of the EOB model, for the DD con- 
figuration. Unlike the nonspinning case the orbital 
frequency fl continues to grow during the plunge. How- 
ever, the EOB light-ring is very close to the peak of the 
numerical /122, as discussed in Sec. Ill Al Note that during 
the ringdown, the frequency computed from the numeri- 
cal /122 shows increasingly large oscillations. We also plot 
the frequency computed from the numerical iff 4 model. 
This frequency shows much smaller, and bounded, oscil- 
lations deep into the ringdown regime. 
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i For Enhanced LIGO, we use the power spectral density given at 
http : //www . ligo . caltech.edu/~rana/NoiseData/S6/DCnoise.txt 
for Advanced LIGO, we use the broadband 
configuration power spectral density given at 
http : //www. ligo . caltech.edu/advLIGD/scripts/ref_des . shtml 



FIG. 7: (color online) We show the amplitude and frequency 
of the numerical and EOB mode /i22, the EOB orbital fre- 
quency and the frequency of the numerical mode for the 
DD configuration. The vertical line marks the peak of the 
amplitude of the numerical waveform. The EOB light-ring 
is 0.3M before the peak and is too close to be shown in the 
figure. 



D. Comparing the gravitational- wave modes hi m 

Here we generate inspiral higher-order modes, hg m , us- 
ing the same dynamics-adjustable parameters calibrated 
to the numerical /122 mode in the previous section. The 
EOB-NQC parameters and the EOB-waveform param- 
eters for these modes are not calibrated, since higher- 
order numerical waveforms show large numerical errors 
before reaching their peaks. For this reason, we constrain 
the comparison between numerical and EOB higher-order 
modes to the inspiral stage. The higher-order modes are 
aligned at low frequencies using the same time and phase 
shifts (modulo a factor of m/2 in the phase shifts) applied 
to the EOB /122 mode. 

In Fig. |H1 we compare the EOB (dashed lines) and nu- 
merical (solid lines) amplitudes of the first five modes 
that dominate the signal power. In the DD case we show 
results only until t = 2000A/ because at later times the 
numerical data are affected by large oscillations, likely 
due to gauge effects. Except for the /132 mode in the DD 
configuration, the agreement is very good for all the sub- 
dominant modes, as well as for the /122 mode, in both the 
UU and DD waveforms. We believe that the difference 
seen for the DD /132 mode is due to the lack of knowl- 
edge of PN spin couplings in the amplitude of the higher 
modes. In fact, only the leading-order PN spin term is 
known in the amplitude of the /132 mode, and no PN spin 
terms are known in the amplitudes of the other subdom- 
inant modes shown in Fig. [5J Being rcsummed in the 
form of Eq. (|30[) , the leading spin term in /132 leaves a 
residual term in fa at the leading order. We test two 
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FIG. 8: (color online). Comparison of the numerical (solid lines), EOB (dashed lines) and Taylor-expanded (dotted lines) 
amplitudes of the dominant and leading subdominant (I, m) modes for the UU (left panel) and DD (right panel) configurations. 
The inset shows the amplitudes for the dominant (2, 2) mode during the late-inspiral and plunge in the DD configuration, 
without the addition of EOB-NQC and EOB-waveform adjustable parameters. 



choices of the odd-parity source term in Eq. ([55]) . L e g 
and H e g, and find this residual term always dominating 
over all the nonspinning terms and causing to de- 
crease and cross zero at high frequency, thus showing the 
odd behavior in the DD (3,2) mode of Fig. [5] We also 
try to apply the p-resummation discussed in Sec. Ill CI on 
the spin terms of the f^2- Although when applying the 
p-resummation, the leading order residual term in ,032 is 
reduced byl/£=l/3 with respect to the residual term in 
^32, it still dominates over other terms and causes (/O32) 3 
to cross zero at high frequency. 

In Fig. [U we also show the Taylor-expanded PN am- 
plitudes (dotted lines) . Their expressions can be read 
from Refs. UU and (24|, and they depend on dynamical 
variables only through the orbital velocity. We calculate 
these amplitudes using the non-Kcplerian orbital velocity 
defined in Eq. (|52"j) for the leading term and the Kcplc- 
rian orbital velocity for all the next-to-leading terms. We 
calculate the non-Keplerian and Keplerian velocities us- 
ing the EOB dynamics. That is to say, these amplitudes 
and the resummed amplitudes are calculated using ex- 
actly the same dynamical evolutions. In particular, the 
energy flux in the dynamics is always modeled by re- 
summed waveforms, even when we calculate the Taylor- 
expanded PN amplitudes. These Taylor-expanded PN 
amplitudes are not to be confused with the amplitudes 
of the adiabatic PN approximants, such as the TaylorTl 
and TaylorT4 approximants (4l| , because the underlying 
dynamics of the latter is completely different. In Fig. [31 
although the Taylor-expanded PN amplitudes work rea- 
sonably well for the /122 mode during inspiral, and prob- 
ably by chance also for the /132 mode in the UU config- 
uration, their performance is not as good as that of the 



resummed amplitudes in general. Especially, for the /144 
and /142 modes, the Taylor amplitudes are not monotonic. 
This unpleasant behavior is caused by their 1PN order 
non-spinning terms. Furthermore, the insert of Fig. [5] 
shows that the performance of Taylor-expanded PN am- 
plitudes becomes worse for the /122 mode during the late 
inspiral and plunge in the DD configuration. 

Given the current information from PN theory and nu- 
merical simulations, we consider the agreement in Fig. [8J 
reasonable and do not dwell further on the choice of the 
waveform modeling options. The differences have little 
impact on the EOB model since the largest difference in 
ft-32 affects the energy flux by less than 10~ 4 , which is 
overwhelmed by other uncertainties in the EOB dynam- 
ics. 

For the five dominant modes, the relative differences 
between the numerical and EOB hg m frequencies are 
within 0.5%, except for the (3, 2) mode in the DD con- 
figuration where the difference is within 1%. Since the 
h£ m frequency depends on both the orbital frequency and 
its amplitude, the larger amplitude difference in the (3, 2) 
mode affects its gravitational- wave frequency. Except for 
the (3, 2) mode in the DD configuration, all frequency 
agreement is within the numerical errors. 



IV. CONCLUSIONS 

In this paper, we carried out the first calibration of 
the spin EOB model to accurate numerical-relativity sim- 
ulations of spinning, non-precessing black-hole binaries. 
We focused on two equal-mass black-hole binaries hav- 
ing spins both aligned, or both anti-aligned with the or- 
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bital angular momentum, and dimcnsionlcss magnitude 
- 0.44 lH. 

For the EOB conservative dynamics, we adopted the 
spin EOB Hamiltonian suggested in Rcfs. [HI, [2(|, aug- 
mented with the 4PN order non-spinning parameter as 
and two adjustable parameters. For the EOB non- 
conservative dynamics, we employed the gravitational- 
wave energy flux which includes spin effects and which 
has been computed using the factorized multipolar wave- 
forms of Ref. 0. 

As in previous cases [IH, UH , we aligned the EOB and 
numerical waveforms at low frequency over a time in- 
terval of 1000M, and minimized the difference between 
numerical and EOB waveforms by calibrating a hand- 
ful of EOB-adjustablc parameters. In particular, in 
this first exploration, we calibrated two EOB-dynamics 
adjustable parameters [b(v) in Eq. (|16bp which intro- 
duces a spin-orbit term at 3.5PN order, and a|g N in 
Eq. (|19p which introduces a 3PN spin-spin term], and 
three EOB-NQC adjustable parameters [see Eq. ([34]) ] 
which enter the gravitational-wave energy flux and the 
EOB gravitational-wave (2,2) mode. Finally, we also cali- 
brated the EOB- waveform adjustable parameter Ai^ atch . 
Quite interestingly, similar to the case of non-spinning 
waveforms, we found that for spinning waveforms, once 
the EOB-dynamics adjustable parameters are calibrated 
at low frequency, the EOB light-ring coincides with the 
peak of the numerical-relativity waveform. Thus, for 
both spinning and non-spinning binary black holes, the 
EOB light-ring marks the most natural point at which 
to match the EOB inspiral-plungc waveform to the EOB 
merger-ringdown waveform. 

In the equal-mass, spin aligned case, we found that 
phase and fractional amplitude differences between the 
numerical and EOB (2, 2) mode can be reduced to 0.01 
rads and 1%, respectively, over the entire inspiral wave- 
forms. In the equal-mass, spin anti-aligned case, these 
differences can be reduced to 0.13 rads and 1% during 
inspiral, and to 0.4 rads and 10% during merger and 
ringdown. The waveform agreement is within numeri- 
cal errors in the spin aligned case while slightly over nu- 
merical errors in the spin anti-aligned case. Despite this 
difference, we found that using Enhanced LIGO and Ad- 
vanced LIGO noise curves, the overlap maximized with 
respect to reference time and phase between the EOB 
and the numerical (2,2) mode, is larger than 0.999 for 
binaries with total mass 30-200M Q . This is well above 
the accuracy requirement of binary black-hole waveforms 
for detection and measurement purposes in gravitational- 
wave observations [40(. 

In addition to comparing the numerical and EOB 
waveforms for the leading (2,2) mode, we also compared 
them for the next four subleading modes. Except for the 
/132 mode in the DD configuration, the amplitude and 
frequency agreements are very good for all the subdom- 
inant modes, as well as for the /122 mode, in both the 
UU and DD waveforms. We believe that the difference 
seen for the DD /132 mode is due to the lack of knowledge 



of PN spin couplings in the amplitude of the subleading 
modes. 

The spin EOB Hamiltonian [lj| [20[ adopted in this pa- 
per was an excellent starting point to explore the calibra- 
tion of the EOB model against spinning numerical sim- 
ulations; however, as discussed above, and in particular 
in Sec. Ill Bl the spin EOB Hamiltonian we used exhibits 
some unusual behaviour. Especially when extended at 
4PN and 5PN, in some regions of the parameter space 
the Hamiltonian does not have an ISCO or the ISCO ra- 
dius grows as the spin magnitude increases. This is op- 
posite to the result in the test-particle limit case. More- 
over, although the spin EOB Hamiltonian has a light 
ring, in some regions of the parameter space (including 
the anti-aligned case discussed in this paper) the orbital 
frequency does not reach a maximum. Those features 
turned out to be crucial when calibrating the EOB model 
to nonspinning numerical waveforms, and we believe will 
be crucial also when modeling spinning numerical wave- 
forms. We found that the lack of those features in the 
current EOB Hamiltonian is due to the ad hoc spin cou- 
pling term H e g pai t(R, P , or) , defined in Eq. (fl~3| . This 
spin coupling term does not reproduce the results of a 
spinning test-particle at PN orders higher than 2.5PN. 
Analyses using an improved spin EOB Hamiltonian [32[ 
obtained by building on Ref. [23[ have shown that those 
features can be recovered. 
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Appendix A: Tortoise coordinate in Cartesian 
implementation 

We start with the definition of the radial tortoise coor- 
dinate given in Eq. ([25)1: dR* / dR = l/£ a (R), The invari- 
ance of the action gives Pr* = P^dR/dR* = Pn£ a (R)- 
In evolving the EOB dynamics, we adopt the dynamical 
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variables R, P* , Si and S^. The transform from P to above, we can write the transform in components as 
P* is a coordinate transform, not a canonical transform. 
In this section, we derive explicitly the transform to tor- 
toise coordinate for the Hamiltonian and Hamilton equa- 
tions of motion implemented in Cartesian coordinates. 

The transform between P and P* is determined by 
the invariance in their tangential components and the 
rescaling in their radial components, that is 
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Choosing three independent equations out of the four or explicitly as 
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In the spin EOB Hamiltonian, we shall replace P with 
T~ 1 P* . The equations of motion for R and P* are 
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where the matrix dP* jdX^ can be written in T and P* 
as dP?/dXi = dT^/dX* (T _1 ) fc 'P;*. 
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